Free-energy Barrier to Melting of Single-chain Polymer Crystallite 
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We report Monte Carlo simulations of the melting of a single-polymer crystallite. We find 
that, unlike most atomic and molecular crystals, such crystallites can be heated appreciably above 
their melting temperature before they transform to the disordered "coil" state. The surface of the 
superheated crystallite is found to be disordered. The thickness of the disordered layer increases 
with superheating. However, the order-disorder transition is not gradual but sudden. Free-energy 
calculations reveal the presence of a large free-energy barrier to melting. 
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It is easy to supercool liquids but difficult to superheat 
solids. The reason is that the surface of a crystal can 
melt at a temperature that is well below the bulk melt- 
ing temperature. Q As a consequence, solids usually melt 
from the surface inward without significant superheating. 
There exist experimental studies that show superheating 
of solids, but in these experiments the crystals are either 
confined in a non-melting matrix Q or the experiments 
reveal superheating of one particular crystal surface only. 

In the present paper, we consider the melting of a macro- 
molecular crystal. The freezing of polymer crystals has 
been the subject of much experimental and theoretical 
work, but much less attention has been paid to the melt- 
ing of polymer crystallites. This is also true for 
the case of a single-chain polymer crystallite. In fact, 
an extensive numerical study of the freezing of such a 
crystal has been reported by Muthukumar and cowork- 
ers |(|fjjand by Fukui et al. M. However, there is no cor- 
responding analysis of the melting of single-chain crys- 
tals, in spite of the fact that this is the simplest ex- 
ample of polymer-crystal melting. The freezing and 
melting of single-polymer crystallites is expected to dif- 
fer from the corresponding phenomenon in systems of 
small molecules. The reason is that the size of a single- 
chain crystallite is limited by the length of the polymer 
chain. This implies that, even if such a crystal would 
have its equilibrium morphology, it would melt well be- 
low the bulk melting temperature. Moreover, when a 
single polymer partially melts (or dissolves), then the 
molten monomers cannot escape from the surroundings 
of the crystallite. Rather they stay around as a "corona" 
and can, in this way, affect the remainder of the melt- 
ing process. The simulations presented below show that 
these features make the melting of polymer crystallites 
qualitatively different from that of atomic or molecular 
crystals. 

To study the melting of a single-chain crystallite, we used 
a lattice model. Lattice models provide a highly simpli- 
fied picture of freezing and melting. Nevertheless, it has 
been shown M that such models are sufficiently flexi- 



ble to account for the phenomenon of surface melting 
in simple "atomic" systems. The polymer lattice model 
that we used in the present study is described in Ref. 
]j~0| |. In this model, polymers live on a simple cubic lat- 
tice, but the polymer bonds can be directed both along 
main axes of the lattice and along the face and body 
diagonals: 26 directions in all. The polymers can be 
semi-flexible and have an attractive nearest-neighbor in- 
teraction. This interaction is anisotropic: parallel poly- 
mer bonds attract more strongly than non-parallel bonds. 
Increasing the anisotropy of the polymer-polymer inter- 
action stabilizes the crystalline state with respect to both 
the molten- globule ("liquid") state and the coil ("vapor") 
state. |l(]] In its most general formulation, this model 
is characterized by three energy parameters. E c is the 
energy required to make a kink in the polymer. The 
larger E c , the stronger the tendency of the polymers to 
be linear. All non- collinear bonds are assumed to have 
the same energetic cost. The interaction between two 
non-bonded adjacent bonds is controlled by two parame- 
ters: the isotropic interaction energy B, which measures 
the energy change when forming one polymer-solvent 
contact from polymer-polymer and solvent-solvent con- 
tacts, and the closely related anisotropic interaction en- 
ergy E p , which only acts between parallel polymer bonds. 
By varying E c , B and E p we can change the "phase- 
diagram" of the polymer. A flexible (E c = 0) polymer 
with a large B undergoes a coil-globule transition. In 
contrast, a polymer with a large value of E p but small B 
will go directly from the coil state to the crystalline state 
- even when E c = 0. In what follows, we consider the sim- 
plest model of a crystallizing polymer, namely one where 
E c and B vanish, but E p does not. At first sight, it would 
seem that such a fully flexible polymer cannot form an 
ordered crystal. However, the interaction energy E p has 
the effect to align polymers in the crystalline state. We 
performed simulations of a single-chain crystallite con- 
sisting of 1024 monomers. To minimize finite-size effects, 
we used periodic boundary conditions and a simulation 
box containing 256 3 lattice points. This was large enough 
to avoid interactions of the polymer with its periodic im- 



1 



age, even in the coil state. 

Initially, we prepared the polymer chain in a rectangular 
crystal with a folding length of 4 monomers and lateral 
dimensions of 16 by 16 lattice units at a temperature 
T = 2.0E p /kB- Subsequently, this crystallite was an- 
nealed for 10 6 MC cycles, where one cycle corresponds 
to one trial move per monomer. This annealing results 
initially in a thickening of the crystallite. Longer anneal- 
ing does not result in further changes of the crystallite 
shape. This suggests that, at the end of the annealing 
run, the crystallite has reached its equilibrium morphol- 
ogy. Having thus prepared an equilibrated single-chain 
crystal, we slowly heat it. To this end, we decrease the 
ratio E p /{k B T) by 0.01 every 10 6 MC cycles. Figure @ 
shows the average potential energy < E/E p > of the sin- 
gle chain as a function of temperature. We choose the 
potential energy of a fully ordered bulk crystal as our 
zero of energy. As can be seen in figure [l], there is a 
sudden change in < E/E p > around 3.125E p /kB- This 
change is due to the sudden and irreversible melting of 
the compact crystallite to an expanded coil state. This 
transition to the coil state is preceded by an apprecia- 
ble pre-transitional increase of the energy. Figure |l| also 
shows the temperature dependence of the fluctuations in 
the internal energy. As is expected, the pretransitional 
rise of the internal energy is accompanied by an increase 
in the energy fluctuations. 




FIG. 1. Stepwise heating curve of the potential energy 
< E/Ep > (spheres) and its fluctuations (triangles) for a 
single-chain crystallite of a lattice polymer with 1024 units. 
The system is allowed to relax during 10 6 MC cycles at each 
temperature. The temperature was increased by decreasing 
E p /(k B T) with a fixed step size of 0.01. Another 10 6 MC cy- 
cles were used to compute average properties at each tempera- 
ture. The dashed line indicates the approximate temperature 
where irreversible melting took place in the simulations. 

In order to distinguish between bonds that belong to 
the crystalline core and those that belong to the corona, 
we need to have an order parameter that is sensitive to 
the environment of a bond. We found that bonds inside 



a crystallite, always have more than 5 parallel neighbors, 
while the overwhelming majority of bonds in the coil 
state have less. We therefore refer to bonds with fewer 
than 5 parallel neighbors as " molten units" . In figure ||, 
we show the dependence of the number of molten units 
on temperature. At low temperatures, about 100 of the 
units that are identified as " molten" , can be attributed 
to the bonds on the hairpin folds at the crystallite sur- 
face. When the crystallite is heated close to the onset 
of abrupt melting, the number of molten units increases 
significantly. 




FIG. 2. Stepwise heating curve of the number of molten 
units (see text) under the same conditions as shown in figure 





To follow what is happening, it is useful to trace the 
structural change on heating. To this end, we analyze 
distribution of local crystalline order as a function of the 
distance from the center of mass. We can quantify the 
degree of local crystalline order by counting the num- 
ber of parallel nearest-neighbor bonds around each bond 
(excluding two adjacent bonds along the chain). In a 
dense, ordered state, the number of parallel neighbors is 
equal to 24. In the coil state, it is more than an order of 
magnitude lower. In figure ^A, we show the local crys- 
talline order as a function of the distance to the center 
of mass of the chain. The figure shows that, as the tem- 
perature increases, the crystalline core is surrounded by 
a "corona" of bonds that occur in a disordered environ- 
ment typical of a polymer coil. Within the time scale of 
10 7 MC cycles, the abrupt melting can be observed at 
about 3.105-Ep/fcs. Figure ||B shows a snapshot of the 
crystallite at a temperature 3.0&6E p /kg, i.e. just below 
the irreversible melting transition. The snapshot indeed 
shows a crystalline core surrounded by a dilute "gas" of 
coil bonds. It should be stressed that these bonds belong 
to many different loops. Hence, the surface disordering 
is not limited to the chain ends. 
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FIG. 3. Radial distributions (3A) of the local crystalline 
order of the bonds (see text) for a series of JcbT/E p values. 
The results are averaged over 5 x 10 4 samples, each with 100 
MC cycles interval. The lines are drawn as a guide to the 
eye. Figure 3B shows a snapshot of single chain at a tem- 
perature 3.086-Ep/fcs, i.e. just below the irreversible melting 
temperature. 



The local crystalline order in the disordered region and 
its range, increase with temperature before the point of 
instability is approached. These results should be asso- 
ciated to the increase of molten units with temperature, 
which enhances the possibility of parallel packing in the 
disordered region as well as its outreach. The uncertainty 
of the tails in the figure ||A is due to the scarcity in statis- 
tics. 

In addition, we found that as long as abrupt melting 
does not take place, the number of molten units can be 
changed reversibly by varying the temperature. Such be- 
havior is expected when there is a free-energy barrier that 
separates the states with partial surface melting from the 
fully disordered coil state. 
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FIG. 4. Free energy barrier for crystallite melting. The 
curves show the dependence of the free energy on the number 
of molten polymer units. The free energies were calculated 
using umbrella sampling with 15 overlapping windows. Since 
our estimation is not for absolute free energy, the curves are 
shifted longitudinally to meet the first minimum at zero line. 

In our simulations, we can directly determine the free- 
energy barrier for irreversible melting. To this end, 
we compute the probability distribution P{n) to find 
n molten units in the system. The Landau free en- 
ergy of a state with this value of n is then given by 
F(n) — —kBTlnP(n). To improve the statistical accu- 
racy of the calculation for those values of n where P(n) 
is small, we used umbrella sampling. [pT| The results of 
this calculation are shown in figure |f|. In this figure, 
we see two minima of the free energy: one corresponds 
to the crystallite with a disordered surface, the other to 
the disordered coil. Figure [| allows us to determine the 
temperature at which coil and crystalline state are in 
equilibrium. At this temperature of about 2.97 E p /kB, 
there is a high free-energy barrier (about 70fcsT) sep- 
arating the crystalline and the molten states. As the 
temperature is increased, the disordered state becomes 
more stable, but the barrier for "explosive" melting re- 
mains much higher than the thermal energy, up to a re- 
duced superheating of 4.5%. Conversely, as a disordered 
coil is cooled down below the coexistence temperature, 
a large free-energy barrier will inhibit spontaneous crys- 
tallization. In this respect, this simple homopolymer re- 
produces some of the aspects of protein folding. If we 
identify the crystalline state with the native state, and 
the coil state with the denatured state, then the present 
model has a pronounced barrier for both folding and un- 
folding. However, we stress that the pathway for melting 
of superheated single-chain crystals is very different from 
the freezing pathway of supercooled chains. 
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